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ABSTRACT 

In previous work stable approximately axisymmetric equilibrium configurations for magnetic 
stars were found by numerical simulation. Here I investigate the conditions under which more 
complex, non-axisymmetric configurations can form. I present numerical simulations of the 
formation of stable equilibria from turbulent initial conditions and demonstrate the existence 
of non-axisymmetric equilibria consisting of twisted flux tubes lying horizontally below the 
surface of the star, meandering around the star in random patterns. Whether such a non- 
axisymmetric equilibrium or a simple axisymmetric equilibrium forms depends on the radial 
profile of the strength of the initial magnetic field. The results could explain observations of 
non-dipolar fields on stars such as the BO. 2 main-sequence star r Sco or the pulsar IE 1207.4- 
5209. The secular evolution of these equilibria due to Ohmic and buoyancy processes is also 
examined. 

Key words: (magnetohydrodynamics) MHD - stars: magnetic fields - stars: chemically pe- 
culiar - stars: neutron - (stars:) white dwarfs 



1 INTRODUCTION 

Convective stars of various types tend to have small-scale, time- 
varying magnetic fields, which is interpreted as dynamo action 
driven by differential rotation and convection/buoyancy instabili- 
ties. Stars which lack convection tend to display large-scale, steady 
magnetic fields. Where we find magnetic fields in upper-main- 
sequence (MS) stars (> 1.5M G ), which are radiative apart from a 
small convective core, they are large-scale and ste ady at lea s t for a s 
lon g as they have been observed (first detection by Ba bcockl Jl947t) . 
see|Mathvs 2001 for a recent review). Many of these stars have 
roughly dipol ar fields, but many have rather more complex fields, 
such as t Sco jDonati et alj2006f) . We find magnetic field s of simi- 
lar geometry and total flux in some white dwarfs (see e.g. lSchmidtl 
1200 ll : IWickramasinghe & Ferrariol 120051: also magnetic fields in 



Type-la SN ejecta - see e.g. Stritzinger & Sollerman 2007). Again, 



many WDs have dipolar fiel ds and many have stro nger quadrupolar 
and octupolar components teeuermann et alj|2007l) , although none 
have very small-scale fields in the same way that convective stars 
do. 

In neutron stars (NSs), there is little o bservational constrain t 
on the precise geometry of the fields (see e.g. lHarding & Lai 2006). 
However, there is evidence in some NSs for a non-dipolar configu- 
ration. For instance, the pulsar IE 12 07.4-5209 has a spi ndown- 
inferred dipole of 2 - 4 x 10 12 G i lPavlov et alj |2002|) while 
ISanwal et al.1 d2002h find a surface field strength of 1.5 x 10 14 G 
from absorption features in the spectrum. Similarly, Bec ker et al.l 
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(2003) find dipole and mean surface field strengths of 10 9 and 10" 
G respectively in the pulsar RBS B 182 1-24. 

None of these non-convective stars display differential rota- 
tion. Given the lack of necessary ingredients for a dynamo, we 
infer that the magnetic fields are 'fossil' remnants in some sta- 
ble equilibrium which therefore evolve not on the relatively short 
Alfven time-scale (years, hours and seconds in MS stars, WDs 
and NSs respectively) but only on longer, diffusive time-scales. 
The fields are stable equilibria which formed from the magnetic 
field left over from the last convective or otherwise 'chaotic' pe- 
riod which the star experienced. In the case of MS stars, this is the 
convective/accretive protostellar phase; WDs (at least, the heavier 
ones which have magnetic fields) evolve mainly from the convec- 
tive cores of intermediate-mass-MS stars; NSs are born out of pre- 
sumably turbulent core-collapse supernovae and experience ~ 100s 
of neutrino-driven c onvection and neutron-finger instability (see 
iBonanno et all2006l) . In addition to their turbulent beginnings, two 
things which all these formerly convective stars have in common 
is the enormous range in field strengths between stars of the same 
class (10 <2 ~ 5 , IO 4 - 9 and 10"" 15 gauss in upper-MS, WDs and NSs 
respectively) and the lack of correlation between the field strength 
and rotation speed (or rather, the inferred rotation speed at birth) 
which is particularly puzzling, given the clear correlation observed 
in active dynamos b etween field strength and rotation speed (see 
e.g. lDibvendull2004l) . 

On the theoretical side, the emphasis has been on finding mag- 
netic field configurations which are in stable equilibrium. Analyti- 
cally, one can produce an equilibrium configuration simply by mak- 
ing sure all forces are balanced and then test its stability to arbitrary 
perturbations either by finding a dispersion relation or, more often, 
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by using an energy method dBernstein et alj[l~958h . Unfortunately, 
it has been easier to demonstrate the instability of various config- 
urations than to find stable configurations. For instance, both ax- 
isymmetric poloidal fields and axisymmetric toroidal fields have 
been shown to be unstable JWrightl l 19731: iMarkev & Tavlerll 19731 ; 
lTavledfl973l : lBrai thwaite 20Q6)). More recently, it has become pos- 
sible to find stable equilibria using numerical methods. Braithwaite 
& No rdlund (2006, hereafter Paper I; see also Braithwaite & Spruit 
|2004) found such a configuration by evolving an arbitrary initial 
magnetic field in time and watching it find its way into a stable 
equilibrium. This equilibrium is roughly axisymmetric, consisting 
of both toroidal and poloidal components (which is logical, given 
that both are unstable on their own) in a twisted-torus configura- 
tion. From the outside, the field looks approximately dipolar, as do 
the fields of many (perhaps most) upper-MS stars and WDs. It is 
the main purpose of this paper to demonstrate the existence - and 
explore the properties and stability - of more complex, non-dipolar 
configurations. 

In Sectional describe numerical magnetohydrodynamic sim- 
ulations where equilibria evolve from arbitrary initial conditions. I 
then use simple analytic methods to explore the stability and shape 
of these equilibria in Section [3] and then look at the slow, quasi- 
static diffusive evolution of these equilibria in Section [4] In Sec- 
tion [5] I discuss the results in the context of astrophysical appli- 
cations and observational consequences before concluding in Sec- 
tion|6] 



2 NUMERICAL SIMULATIONS 

In Paper I it was found that an arbitrary initial magnetic field can 
evolve into a roughly axisymmetric twisted torus configuration. 
Here, the aim is to investigate whether that is the only stable equi- 
librium available or whether there are other equilibria which can be 
reached from different initial conditions. 



2.1 The numerical model and setup 

The setup of the simulations was described in some detail in Pa- 
per I; I give a brief summary here. The star is modelled as a ball 
of self-gravitating ideal gas (ratio of specific heats y = 5/3) of 
mass M, arranged in a polytrope of index n = 3, so that specific 
entropy increases with radius and the star is stably stratified. [This 
has important consequences; the behaviour of a magnetic field in an 
isentropic star will be explored in a forthcoming paper.] The star is 
contained in a computational box of side 4.2R, where R is the stellar 
radius. Surrounding the star is an atmosphere of low electrical con- 
ductivity, which behaves like a vacuum in that the magnetic field 
there relaxes to a potential (current-free) field. The initial magnetic 
field is chosen to resemble that expected to be present at the end of 
a period of convection. In other words, we want to create a chaotic, 
small-scale field. To do this, the initial magnetic field is calculated 
from a vector potential (thus ensuring V • B = 0) which contains 
wavenumbers with a flat power spectrum up to a magnitude k mdx . 
The vector potential is tapered so that 

B~p", (1) 

where p is the gas density and p is a free parameter. Note that if 
the star forms out of a uniformly magnetised cloud and the same 
fraction of flux is lost from all fluid elements, then we expect p = 
2/3. In all cases, some magnetic flux goes through the surface of 



the star; a greater amount for smaller p (see Figure [8] where the 
solid lines represent initial conditions). 

The code used is the stagger code (Nordlund & Galsgaard 
1 19951: iGudiksen & Nordlund 120051) . a high-order finite-difference 
Cartesian MHD code which uses a 'hyper-diffusion' scheme, a 
system whereby diffusivities are scaled with the length scales 
present so that badly resolved structure near the Nyquist spatial 
frequency is damped whilst preserving well-resolved structure on 
longer length scales. This, and the high-order spatial interpolation 
and derivatives (sixth order) and time-stepping (third order) in- 
crease efficiency by giving a low effective diffusivity at modest res- 
olution (128 3 here). The code includes Ohmic and well as thermal 
and kinetic diffusion. Using Cartesian coordinates avoids problems 
with singularities and simplifies the boundary conditions: periodic 
boundaries are used here. 

As always with numerical work, one has to ensure that the rel- 
evant time-scales, which are separated in reality by many orders of 
magnitude, are also sufficiently separated in simulations. Here, the 
three relevant time-scales are the sound-crossing time, the Alfven 
time-scale and the diffusion time-scale, with the ordering 

t s « t a « r d . (2) 

The short sound crossing time ensures that the star remains in pres- 
sure equilibrium, or rather, that it evolves on an Alfven time-scale 
in quasi-static pressure equilibrium, where I define the Alfven time- 
scale as t a = R s/M/2E where E is the magnetic energy in the star. 
Likewise, once the magnetic forces are balanced, the star evolves 
on the diffusive time-scale in quasi-static MHD (and pressure) equi- 
librium. The strength of the magnetic field in the simulations must 
therefore be low enough that the Alfven speed is everywhere much 
less than the sound speed: the ratio of thermal to initial magnetic 
energy is thus chosen to be 400. However, the magnetic energy 
is expected to fall significantly, and to stop the Alfven time-scale 
from becoming comparable to the diffusion time-scale, the strength 
of the magnetic field is artificially boosted as it decays so that the 
total magnetic energy stays constant. To do this, a term is added 
to the induction equation dB/dt = .... + «B, where a is calculated 
at every timestep. A record is kept of the degree of boosting so 
that the 'actual' magnetic energy can be retrieved - it is this quan- 
tity E' which is plotted in Figure [7] where E' is calculated from 
dE' /dt = -2aE' . This holds E and therefore ta constant and main- 
tains the time-scale ordering in l(2j. 

2.2 Formation of equilibria 

The evolution in time of the magnetic field is followed in the sim- 
ulations. The value of & max was set to 36i? 1 and various values of 
p were tried: 0, 1/3, 2/3 and 1. In all cases, a stable equilibrium is 
reached after a few Alfven crossing times, but both the geometry 
and energy of the equilibrium is found to vary. Two types of equi- 
librium are found: simple axisymmetric (as found in Paper I), and 
more complex non-axisymmetric. 

Figure Q] shows the initial conditions (of the p = 1 run). Note 
that the star takes up only a small part of the box - this wastes some 
computational power but is necessary to ensure an approximately 
correct treatment of the potential field outside the star. For the run 
with p = 1, the magnetic field after around 9 Alfven crossing times 
is displayed in Figure|2] The configuration is roughly axisymmetric 
and contains both toroidal and poloidal components. 

When p = 0, a non-axisymmetric equilibrium is reached. This 
is displayed in Figure [3j the figure shows the star from both sides. 
The configuration consists of twisted flux tubes below the surface 
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Figure 1. The initial conditions: the top left shows the whole computational 
domain; the main picture shows just the star, which is represented by blue 
shading. 




Figure 2. The simplest of the equilibria, which forms when the initial mag- 
netic field is sufficiently concentrated into the centre of the star. The config- 
uration is roughly axisymmetric, with a twisted flux tube looped around the 
equator and a poloidal field passing through the middle of the star. 



of the star. The flux tubes lie horizontally and meander around the 
star at some distance below its surface. 




Figure 3. A non-axisymmetric equilibrium viewed from opposite sides of 
the star. The equilibria consist of one or more twisted flux tubes buried just 
below the surface of the star. This kind of equilibrium forms if the initial 
magnetic field is not concentrated into the centre of the star. In making this 
figure, the 'opacity' of the star (blue shading) was set so that the field lines 
just below the surface can be seen but that lines further inside cannot, as a 
visual aid. 



3 STABILITY OF NON-AXISYMMETRIC FIELDS 

In this section, the structure and stability properties of the non- 
axisymmetric field will be examined, and predictions made about 
how the nature of the equilibrium depends on the conserved param- 
eters energy, helicity and flux. 
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Figure 4. Cross-sections of a twisted flux tube below the surface of a star 
of radius R. Above, the axisymmetric case where the flux tube is a circle 
around the equator; below, the non-axisymmetric case where the flux tube 
meanders around the star in some more complex fashion. The directions of 
B\ and B r are indicated by the arrows in the lower part of the figure; B t is 
directed into the page. The poloidal field lines are marked with arrows; it 
can be seen that the boundaries of the wedge are magnetic surfaces. The 
toroidal field (direction into/out of the paper, shaded area) is confined to the 
poloidal lines which are closed within the star. 



3.1 The structure of a twisted flux tube 

Consider a twisted flux tube of length s and width (at the surface 
of the star) aR whose axis lies at some constant radius inside a star 
of radius R, so that a is the angle subtended by the flux tube at the 
centre of the star. The cross-section of such a flux tube is illustrated 
in the lower part of Figure [4] It is assumed that a <k 1 (see sec- 
tion [37IJ. The three components of the magnetic field are toroidal 
(parallel to the length of the flux tube), radial and latitudinal; av- 
erage values are denoted respectively by B t , B\ and B r . The radial 
and latitudinal components will be collectively referred to as the 
poloidal component. [Note that this coordinate system is local to 
each position along the flux tube.] Since the flux tube cannot have 
any ends, it must be joined in a loop. This loop can be wrapped 
around the star in any fashion; it is assumed below that the radius 
of curvature is significantly greater than the tube's width. 

We can immediately compare this picture to the simulations. 
In Figure[5]we have five flux-tube cross-sections from the run with 
p = at time t = 17t a . The upper three frames are centred on 
different parts of one flux tube; the lower two frames on different 
flux tubes. In each case, the axis of the flux tube is perpendicular 
to the page. Toroidal field is represented by blue/red shading and 
poloidal field is represented by arrows as well as by contours of the 
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Figure 5. Cross-sections of flux tubes in the run with p = at time 
t = 17ta. The curved grey line towards the right is the surface of the star 
and the centre of the star is on the left. The blue/red shading represents the 
toroidal field component (out of/into the page) multiplied by m, the horizon- 
tal coordinate in these cross-sections, the cylindrical radius. The poloidal 
component is represented by the arrows and by contours of its scalar poten- 
tial. It can be seen that the arrows are very nearly parallel to the contours 
of the scalar potential, showing that the length scale of variation in the di- 
rection perpendicular to the page is much greater than the length scale in 
the plane of the page, i.e. that the flux tubes meander around the star over 
scales much greater than their width. 
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scalar potential if/ (see e.g. equation 3.6 of lMarkev &TavleJll973h 
defined as 



— = -mB m and — = wB- 

oz am 



(3) 



Here, it is more correct to use cylindrical coordinates w and z, as 
opposed to the spherical r and 1 above, although the two sets are the 
same in the limit of small a. The axis of the flux tube lies on (and 
in fact defines the location of) both the z and 1 axes, so that both 
these coordinate systems are local, as opposed to the usual cylin- 
drical system for instance where z has the same direction at every 
location. Now, this scalar potential is only meaningful if the diver- 
gence of the poloidal field is negligible, i.e. that the relevant length 
scale along the flux tube (i.e. its radius of curvature) is much greater 
than the tube's width. That the arrows in Figure[5](representing the 
real direction of the poloidal field) and contours of if/ are roughly 
parallel confirms that this is a reasonable app roximation. 

Mestel (1961, see also lRoxburgnH r 966) found that in a star 
with an axisymmetric field, the condition that the azimuthal com- 
ponent of the Lorentz force vanishes (necessary since it cannot be 
balanced by the pressure gradient) gives the condition 



B p • V(crfl t ) 



■ 0. 



(4) 



where B p is the poloidal field, B t is the azimuthal (toroidal) com- 
ponent and bt is the cylindrical radius. In other words, poloidal 
field lines are contours of crS t . In the non-axisymmetric case the 
same argument applies: ru in this context is the distance from the 
line passing through the centre of the star, perpendicular both to 
the flux tube axis and to a line connecting the centre of the flux 
tube to the centre of the star (i.e. the dotted line on the left of 
Figure [4}. It is also shown in Figure [4] that the toroidal compo- 
nent of the field is confined to the volume enclosed by the largest 
poloidal line closed within the star, and is zero outside this volume. 
The reason for this is that long-lasting currents cannot exist outside 
the star and consequently the toroidal field outside the star has to 
be zero, and that tuBj, has to have the same value along a given 
poloidal line. A hand-waving argument as to how the toroidal field 
outside the largest closed poloidal line is destroyed is as follows. A 
poloidal line which crosses the stellar surface would 'unwind' if it 
contained any toroidal component. In this process, toroidal field is 
redistributed along the whole of the poloidal line. Toroidal field in 
the atmosphere is then destroyed by reconnection on a short time- 
scale, and the toroidal field inside the star is redistributed along the 
poloidal line to where it has been lost. This occurs on the Alfven 
time-scale. In this way, all toroidal field on the poloidal line is 
quickly transferred into and destroyed in the atmosphere. By anal- 
ogy, one can think of twisted elastic bands: those which are con- 
nected up in a loop stay wound, while those with free ends unwind. 

These two properties are consistent with the results from the 
simulations described in the previous section - in Figure [5] we see 
that the poloidal field lines do roughly coincide with contours of 
mB t , and that the toroidal field is very weak outside of the last 
closed poloidal line. In the figure, the contours are evenly (linearly) 
spaced, and weak B t of either direction is coloured white. 



3.2 Reaching equilibrium: energy considerations 

The volume of the flux tube is saR 2 /3 and the magnetic energy it 
contains is 



f B 2 



saR 2 1 
3 Sn 



Bf + Bt + 1 + 



-> 2 



B. 



(5) 



where s is the length of the tube (the length on the surface of the 
star), a is the angle subtended at the centre of the star and R is the 
radius of the star. The first two terms in the square brackets and 
the first term in curved brackets represent the field energy inside 
the star from each of the three components, and the last term rep- 
resents the field energy outside the star. The field outside the star is 
potential (current-free) so can be calculated by solving the Laplace 
equation, given B, at the surface. The form of B, at the surface is 
not known, so I have simply assumed a basic sinusoidal form with 
the coefficient a, (of order, but less than unity) relating the mean 
radial field component inside the flux tube to the radial field com- 
ponent on the surface. It will be seen below that the energy above 
the surface is not important unless a is very small. 

We now allow the the length and width of the flux tube to 
change adiabatically. The flux tube will stretch or contract until 
dE/da = 0, i.e. an equilibrium is reachecfl Now, we are consid- 
ering a star where the magnetic pressure is everywhere very much 
less than the thermal, so magnetic forces can have only a negligi- 
ble effect on the gas pressure. Also, the star is stably stratified so 
that energy is required to move fluid elements in the radial direction 
since work has to be done against the buoyant restoring force, and 
again, the magnetic forces are weak in comparison to the buoyancy 
force. We therefore have the following restrictions on the displace- 
ments £ which take place during this adjustment to equilibrium: 



V • f =s and r • f * 0. 



(6) 



In other words, weak magnetic forces have to avoid doing work 
against stronger forces. These approxi mations are standard in anal- 
ysis o f MHD instabilities in stars (e.g. lPitts & Tavledfl986l : ISpruitl 
2002). The consequence of these restrictions is that as long as only 
reversible processes are allowed, the surface-field factor a r stays 
constant as does the volume of the flux tube, since no matter can 
leave or enter it, which can be expressed by stating that the area of 
the flux tube on the stellar surface F = saR = const. The radial 
distribution of magnetic energy also cannot change since gas and 
therefore the magnetic energy frozen into it cannot be transported 
in the radial direction. In addition, we have: 



ainfi, 



-1 



dlnfl. 



1 and — = 0, 



(7) 



d In a din a din a 

which follow simply from flux-freezing as well as from {6]l. We can 
now differentiate ((5} with respect to a and set dE/da = 0, giving 



B:*Bt + 



i 2 

Aa-a , 
An 



-B »*H 1 + 



4na 



(8) 



where we have noted from the geometry that B\ as aB r to arrive at 
the part on the right hand side. Remembering that a r is somewhat 
less than 1, we now look at the case where a is large enough so that 
the term in <[8j containing it can be ignored. We are left with 



B t ~ B[ ~ aB r . 



(9) 



Therefore the adjustment to equilibrium consists in stretching 
or contracting until the toroidal and latitudinal components are 
roughly equal. This principle is also (partially) applicable to ax- 
isymmetric equilibria - see Section |572l 

On the other hand, if a is very small, l[9) becomes instead 

4na 

Bi * aB r . (10) 



3a 2 



1 It is straightforward to verify that d 2 E/da 2 > at this point, i.e. that the 
equilibrium is stable to stretching/contracting perturbations. 
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However, at least in Ap stars there is no observational evidence for 
this, i.e. for magnetic fields with very small-scale structure. In any 
case, a small-scale field would evolve and decay on a timescale 
shorter than the stellar lifetime. (|9j will be used below instead of 
dipt , although it should be noted that this point might conceivably 
need to be revisited in future. 



3.3 Magnetic flux and helicity 

If we define a vector potential A from B = V x A, it can be shown 
dWoltieJ[l958h that the magnetic helicity H = JB • AdV, which 
is a global rather than local quantity, is conserved in the limit of 
infinite conductivity. 

Consideration of the length scales R, aR and R in the toroidal, 
latitudinal and radial directions respectively and from the definition 
of A we see that 

g r -- — , S, ~ — and B r -, (11) 

aR R R ' aR 

which using (|9) becomes A, ~ A\ ~ A r /a ~ RB t . The helicity of 
the flux tube is now given by 



axis of symmetry 

A 



AdV ~ 2FB r O t « 40„cD t 



(12) 



where <D t « (l/2)aR 2 B t and O p » (l/2)saRB r = (l/2)FB r are the 
toroidal and poloidal fluxes; note that all coefficients of order unity 
are rather approximate. The helicity, as well as both O t and O p , 
are conserved during adiabatic changes in a and s and evolve only 
due to diffusive processes; this is described in Section [4] It will 
be seen in Fig [7j that helicity is indeed roughly conserved during 
the formation of equilibria on the Alfven timescale, so remains a 
useful tool to aid our understanding of dynamical relaxation into 
equilibrium. 

Using <j9j, we can now give a as a function of the conserved 
quantities defined above: 

2 2 <D t F O t 



(13) 



so that if we know the area F and the ratio of toroidal to poloidal 
fluxes of a tube, we can calculate its equilibrium width aR and 
length s. 

The above arguments can be applied not only to the flux tube 
as a whole, but also to any small length 6s of the tube with area 
SF. Along the length of the flux tube O t must be constant, but B r a 
26<b p /6F can vary. There is no reason why a single flux tube needs 
to have a uniform width - at every point along its length, there is 
a local equilibrium giving a as a function of the toroidal flux and 
the local value of B r (remember from (|7) that B, is constant during 
adiabatic changes in a). This can be seen in the first three plates in 
Figure|5] which are centred on the same flux tube at different points 
along its length - at the position of the third plate, a is greater than 
at the other two. 

Neighbouring flux tubes need to have their poloidal compo- 
nents in opposite senses in order to avoid a discontinuity at the 
boundary between them. If the neighbouring tubes are actually 
parts of the same tube, then the toroidal components will clearly 
also be in opposite senses, i.e. in the cross-section, if the toroidal 

2 Magnetic helicity is gauge-independent only if the boundaries of the re- 
gion of integration are either periodic or magnetic surfaces with B ■ dS = 0. 
The boundaries of the flux tubes here (Figure f4) are indeed magnetic sur- 
faces; note that B goes to zero at infinity. 




Figure 6. The cross-section of a star containing the simplest imaginable 
zero-helicity equilibrium. The shaded areas represent the toroidal compo- 
nent of the field - dark (light) shading represents field going into (out of) 
the page. The vertical axis is the star's axis of symmetry. 



field of the tube goes into the page, the toroidal field of the neigh- 
bouring tube section will be directed out of the page. If the neigh- 
bours are part of other tubes, the poloidal components still need 
to have opposite senses but there is no constraint on the direction 
of the toroidal part. In fact, there is no way that one flux tube can 
'know' about the direction of the toroidal field in any of the other 
tubes present, because they are insulated from one another by re- 
gions of zero toroidal field. In the first plate of Figure[5]we see that 
the two tube sections above and below the tube in the centre have 
both poloidal and toroidal fields in opposite senses to the tube in the 
centre, so that all three have helicity (which is essentially just the 
product of the toroidal and poloidal fluxes) of the same sign, and 
may in fact be different parts of the same tube. In the second plate, 
the lower tube has opposite helicity to the other two tube sections 
visible. This means then that the equilibrium has lower net helicity 
than an otherwise identical configuration with the toroidal field of 
this tube reversed; if we also reduced the strength of this new con- 
figuration we could make a field of equal helicity but lower energy. 

Therefore, there are many stable equilibria with the same he- 
licity but different energies, and it is possible that an arbitrary ini- 
tial field will evolve on the Alfven timescale, whilst conserving 
helicity, into a higher-energy, local equilibrium rather than into the 
lowest energy state for its helicity. Stable configurations with zero 
net helicity are possible in principle. The cross-section of a simple 
zero-helicity field is illustrated in Figure(6] Whether such a higher- 
energy local equilibrium can actually be reached from realistic ini- 
tial conditions is another question. Also, it is plausible that diffusive 
evolution (where helicity is no longer conserved) on a much longer 
timescale may coax the configuration into a lower-energy equilib- 
rium. 



3.4 Simulations: a quantitative analysis 

It is now informative to have a more quantitative look at the simu- 
lations. First, one can measure the magnetic helicity and magnetic 
energy. In Figure [7] these are plotted for the runs with p = 0, 1/3, 
2/3 and 1 (p was defined in (Q3). Time, on the horizontal axis, is 
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5 10 15 20 

t / r A 

Figure 7. Log magnetic energy E and helicity H/8rrR (thick and thin 
lines respectively) against time, for runs with various values of the central- 
concentration parameter p. Helicity has been divided by 8ttR to give it units 
of energy and so that the ratio of the two is equal to the helicity length as a 
fraction of the stellar radius, ru/R. 

given in terms of Alfven times ta = R -\JM/2E where E is the mag- 
netic energy. All runs have the following in common: the helicity 
falls throughout at some low rate, while the energy falls signifi- 
cantly during the first couple of Alfven time-scales before settling 
down to a similar rate of decay as the helicity. This can be inter- 
preted as the formation of a stable equilibrium, followed by decay 
of the equilibrium configuration on a longer, diffusive time-scale. 
Differences between the runs are also visible in the figure: once 
equilibrium has been reached, the runs with p = 1 and 2/3 experi- 
ence a slower drop in energy and helicity than the runs with lower 
values of p. This is presumably because the length scales in the ax- 
isymmetric configurations which form in the p = 1 and 2/3 runs 
are greater and the diffusive time-scales correspondingly greater. 

We can also look at the radial profile of the magnetic energy 
density B 2 /Sn to see where most of the energy resides. In Figure[8] 
are plotted radial profiles of the magnetic energy density at differ- 
ent points in time for each of the four runs with different values of 
p. We can see that the energy profiles change little over the Alfven 
time-scale when the equilibrium is forming, but that on a longer 
time-scale it becomes shallower and eventually the energy density 
peaks at around r » 0.8R, where the axes of the flux tubes are lo- 
cated. Transport of flux in the radial direction is hindered by the sta- 
ble stratification and can only take place on a diffusive time-scale. 
Precisely which kind of diffusive process (magnetic or thermal) is 
responsible is discussed in Section|4] 

It is also useful to measure the length of the flux tubes. The 
easiest method is to measure the length of the lines on the surface 
where B r - 0, as the flux-tube axes always lie directly below these 
lines. Figure [9] is a plot of W, the total of the lengths of all flux 
tubes in terms of the stellar circumference, against time, for the 
runs with p = 0, 1/3, 2/3 and 1. Clearly, the runs with higher p 
values, which form axisymmetric equilibria, should have values of 
« 1 because the flux tube is essentially a great circle around the 
The non-axisymmetric equilibria have higher values, because 
the flux tubes are longer and meander around the star in a more 

3 In fact, what I refer to here as 'axisymmetric' is only approximately so. 
The equilibria are generally slightly offset from the centre of the star and 
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Figure 8. Radial profiles of magnetic energy density B j%n (log, as a frac- 
tion of the mean) at various times for the four runs with different values 
of p. In each plot the lines correspond to the times (in Alfven time-scales) 
indicated in the boxes. It can be seen that the profiles change little during 
the formation of the equilibria, and change only on the longer, diffusive 
time-scale. 
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Figure 9. The length of the B, = lines on the surface of the star, divided 
by the circumference of the star (denoted by W), against time, for runs with 
various values of the central-concentration parameter p: 0, 1/3, 2/3 and 1, 
For the higher values (centrally-concentrated initial field) a dipolar field is 
formed, while for the lower values of p a more complex equilibrium results. 

complex fashion; amongst these equilibria, a lower value of p leads 
to longer, narrower flux tubes. If we assume that the whole volume 
of the star is taken up with flux tubes of similar width then a x 2/W 
in the limit of small a. 

The energy in the radial component of the field as a fraction 
of the total energy, E r /E, is plotted in Figure[l0] In light of l(9]l we 

may be a little warped in shape. This is presumably a result of the random 
nature of the initial conditions. 
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Figure 10. The energy in the radial component of the field B, as a frac- 
tion of the total, against time for runs with various values of the central- 
concentration parameter p: 0, 1/3, 2/3 and 1. As expected, the non- 
axisymmetric equilibria have greater radial energy fractions. 



expect that E r /E « 1/(1 + 2a 2 ), although {9]l is really only valid if 
the radial energy profile is fairly flat. If we look for instance at the 
p = run at 1 = 17t a (when the cross-sections in Figure [5] were 
taken), we see that from Figure [9] that W ~ 4.5, so that a « 0.44 
and we expect E r /E ~ 0.72, which corresponds approximately to 
the value of as 0.75 in Figure [TOl 



4 DIFFUSIVE EVOLUTION OF EQUILIBRIA 

First, let us assume that the electrical conductivity is high enough 
and that thermal conductivity is low enough that the relevant dif- 
fusive time-scales are much larger than the Alfven time-scale, i.e. 
that (aR) 2 /n » aR/v A and (crR) 2 /k » aR/v A , where n and k are 
the magnetic and thermal diffusivities respectively and va is the 
Alfven speed. Then the flux tube remains in quasi-static equilib- 
rium described by (9j and d 1 3I > while it is evolving diffusively. This 
is the case in the simulations, as we can see from Figure[7]that the 
equilibrium evolves diffusively over many Alfven time-scales. Fur- 
thermore, assume that the area F of the flux tube remains constant. 
This is justified if the flux tube(s) have filled up all of the available 
space on the surface of the star, i.e. that F = 4nR 2 ; the simula- 
tions show this to be roughly correct, although there are some small 
'empty' spaces, one of which is visible in the last plate of Figure[5] 

In the simulations, we see that the total magnetic energy in 
the equilibria falls, and that the energy moves to larger radii (Fig- 
ures|7]and[8]l, both taking place on some diffusive time-scale (rather 
than on the Alfven time-scale on which the equilibrium forms). 
The overall fall in energy (and helicity, the product of toroidal and 
poloidal fluxes) can be understood as an obvious consequence of 
finite conductivity, which can also explain the flattening of the ra- 
dial energy profile; what is more difficult to explain is the evolution 
from a flat profile into a state where the energy density is higher at 
larger radii than in the centre of the star. Nor can finite conductiv- 
ity easily explain why a should fall - in Figure [9] we see that the 
length of the flux tubes is gradually increasing, and a is therefore 
falling. From d!3t we can infer from a falling a that the toroidal 
flux O t is falling faster than the poloidal flux O p . These phenomena 
can perhaps be better understood in terms of thermal diffusion. 

Nevertheless let us first look at the effect of Ohmic dissipa- 



tion on the flux tubes. Finite conductivity causes toroidal field to 
leak out of the volume enclosed by the largest closed poloidal field 
line, onto poloidal lines which go through the surface of the star. 
Via the mechanism described in Section [3Tl this toroidal field is 
quickly destroyed. The effect of finite conductivity on the poloidal 
component is that gradients in the radial and latitudinal directions 
are smeared, causing the overall amplitude to go down. Also, since 
the direction of the poloidal component is changing over shorter 
length scales (in the horizontal direction) deeper in the star, this 
could explain why the field strength deeper inside the star drops 
faster than that further out. Also, perhaps since the toroidal field is 
confined to a smaller space than the poloidal, it should decay faster 
since the diffusion time-scale is shorter for smaller length scales, 
thus explaining the fall in time of a. 



4.1 Buoyancy 

A more convincing explanation for the rise of flux tubes to the sur- 
face and for the greater loss of toroidal than poloidal flux involves 
thermal rather than magnetic diffusion. [There is an analogous pro- 
cess, in the case where the stratification is due not to an entropy 
gradient but to a composition gradient and the relevant diffusivity 
is that of chemical elements rather than heat. Here, I shall illustrate 
the principle only for the case of entropy stratification and ther- 
mal diffusion.] This process can cause the field to rise towards the 
surface of a stably stratified star in the following way. Imagine a 
magnetised region surrounded by non-magnetised plasma. It must 
be in pressure equilibrium with the surroundings, so that the sum of 
its magnetic and thermal pressures is equal to the thermal pressure 
outside. If it has the same temperature as the surroundings, it will 
be less dense, and rise. In a stability stratified star (specific entropy 
increasing with radius) the magnetised region will cool adiabati- 
cally as it rises, its temperature falling faster than the ambient tem- 
perature, and it will eventually reach an equilibrium where its total 
pressure and density are equal to those of its new surroundings. 
However, it will now have a lower temperature than the surround- 
ings, and any thermal diffusion will cause it to absorb heat and rise 
further into a new equilibrium. Now, this mechanism should cause 
our flux tubes to rise and even in the absence of magnetic diffu- 
sion inside the star, could lead to destruction of the toroidal flux 
as poloidal lines breach the surface of the star and their toroidal 
component is lost to reconnection processes in the atmosphere. 

Consider a flux tube of radius a, containing toroidal (axial) 
and poloidal field components B t and B p surrounded by an unmag- 
netised gas of pressure P and temperature T. The tube is thin, so 
a <K Hp where H P is the pressure scale height. The pressure exerted 
by the field on the surrounding gas, which is equal to the difference 
in external and internal thermal pressures, is given by 



i a i ff + jg 

AP = net - 

Tjta da I %ji 



8;r 



(14) 



where the expression in brackets represents the magnetic energy 
per unit length of the tube. Note that the poloidal component does 
not contribute, which comes from the fact that adiabatic expansion 
of the tube does not change the energy of that component. Denoting 
quantities internal and external to the tube by i and e, we have at all 
times 



and 



Pi + AP. 



Assuming an ideal gas equation of state gives us AT /T 
and the heat flow into the tube per unit length is 



(15) 
AP/P, 
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(16) 

where = R//-i is the gas constant. After some algebra, the vertical 
speed v of the flux tube is found to be 



H P T " sc 



2k 

n 2R. 



v., 



V a - V 



where fl, = SnP/Bf » 1 and 



V = 



dlnT 
dlnP 



and V a = 



dlnT 
dlnP 



(17) 



(18) 



are derivatives of the surroundings and of the equation of state at 
constant entropy, respectively. Note that the speed becomes infi- 
nite as the temperature gradient approaches the adiabatic gradient 
where convection sets in; in reality the speed would be limited by 
factors not considered here such as aerodynamic forces. Except 
near the boundaries of a stably-stratified zone, we can treat the ex- 
pression in brackets in d!7t as a factor of order unity. 

Strictly speaking, the flux tubes in these non-axisymmetric 
magnetic equilibria are not surrounded by a non-magnetised 
medium, although there are gaps between the toroidal-field regions 
of neighbouring flux tubes. Nor is a <k H P \ however, it seems likely 
that expression J17b is still applicable here, albeit only rather ap- 
proximately. More quantitative, local analysis and simulations (per- 
haps also in 2-D) could shed light on this process and on the speed 
at which the flux tubes rise. Very approximately, we have 



T nsc ~ (aRfpiK ~ CTPTkh 



(19) 



where /3 is the ratio of thermal to magnetic energy and tkh is the 
Kelvin-Helmholtz (thermal) time-scale. Whether this is longer or 
shorter than the Ohmic diffusion time-scale depends on the type 
of star under consideration (see Section I5.lt : in the simulations 
presented here, thermal and magnetic diffusivities are equal and /? 
ranges from ~ 10 4 in the centre of the star to ~ 30 near the surface. 



4.2 Horizontal movements 

As the length and width of the flux tubes adjust quasi-statically 
to the loss of toroidal flux, their arrangement under the surface of 
the star can be expected to change, with tubes pushing each other 
around in the horizontal direction. These movements will be slow 
for the most part, but fast readjustments may also occur, particu- 
larly if stress builds up between tubes and a different arrangement 
of flux tubes becomes energetically accessible. The observational 
consequences of this are not very clear; in addition to a change in 
the geometry of the magnetic field one might expect some release 
of energy accompanied by an apparent perturbation to the rotational 
phase (and also to the rotation period, if the moment of inertia is 
changed). 



4.3 The depth and cross-sectional shape of the flux tubes 

From J X 3 1 > we can predict the width of a flux tube every point along 
its length if we know its toroidal flux and the local value of B,. 
However, we still have no argument to explain the cross-sectional 
shape of the tube, or rather, the shape of the part where the toroidal 
field resides. Even in the same simulation at the same time-step, we 
see a variety of shapes, ranging from circular to very elongated in 
the radial direction (see Figure[3), whilst the depth of the tubes does 
appear fairly uniform. It looks therefore that the radial distribution 
of toroidal flux, and very likely the radial distribution of energy, are 



globally fairly uniform and that only the flux-tube width a varies 
from place to place. This uniformity could be a result of the stable 
stratification and the restrictive effect it has on the movement of 
fluid and therefore magnetic flux in the radial direction, in contrast 
to movements in the horizontal direction caused by changes in the 
width Ra of the flux tubes. We therefore expect the radial distribu- 
tion of energy in the newly-formed equilibria to be similar to that 
in the initial conditions. 

Once the initial equilibrium has formed, the configuration con- 
tinues to evolve on a slower, diffusive time-scale as discussed in the 
previous section, and as this happens, the radial distribution of en- 
ergy changes (Figure[8j. After some time, when the flux tubes have 
risen sufficiently and the magnetic energy density is low in the cen- 
tral part of the star, we can define a depth fR above which most of 
the energy is contained, and redefine B t , B\ and B r as the average 
values in that volume within fR of the stellar surface. (|9) becomes 



B t at B[ a yB r , 



(20) 



since the energy argument in Section 13.21 is unchanged, but the 
different geometry changes the relation between the two poloidal 
components. Equations d 1 1 b -J 13b are changed only by factors of 
order unity. If we now let / change while keeping a and s fixed, we 
find that 



a in a 



-l , 



ains 



-1 and 



dlnfi,. 







din/ din/ din/ 

and that the equilibrium condition dE/df = yields 
fl? * B\ - B\. 

Combining with ((20} we have / as y/2a and 

B t * B\ a V2B r . 



(21) 



(22) 



(23) 



Clearly the depth / of the flux tube cannot change adiabatically on 
an Alfven time-scale as a can, but one at least expects that diffu- 
sive processes will tend to produce flux tubes of roughly circular or 
somewhat radially-elongated cross-section. This process can also 
be understood in another way. In general, a circular flux tube can 
be force-free if the hoop stress of the poloidal field is balanced 
by the force resulting from the gradient in the toroidal field. If the 
tube is now squashed into a tube of elliptical cross-section whilst 
maintaining its cross-sectional area, the hoop stress of the poloidal 
component is still equal on the major and minor axes. However, the 
gradient of toroidal field is stronger now along the minor axis than 
along the major, resulting in a restoring force pushing the tube back 
into its original shape. A flux tube can only be force-free if it has 
circular cross-section, corresponding to the lowest energy state for 
a given helicity. 

Looking at the simulations (Figure [5} we see that most of 
the flux tubes do indeed appear to have / similar to or somewhat 
greater than a. However the fourth frame of the figure has / signifi- 
cantly greater than a. This probably comes from effects not consid- 
ered here such as pressure from neighbouring flux tubes or details 
of the buoyancy mechanism which brings the tubes upwards. 



4.4 Axisymmetric equilibria 

In Paper I we found that an axisymmetric equilibrium gradually 
moves outwards, the torus expanding. The rate at which the field 
on the surface increases was found - in the case where the field 
is deeply buried inside the star - to be proportional to the mag- 
netic diffusivity, which works simply to smooth out any gradients 
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in field strength so that any centrally-concentrated field becomes 
more spread out, and therefore become stronger nearer the surface 
although the total energy of the field is falling 4 ]. However, in Paper I 
we did not look at the effect of the two types of diffusivity when the 
axis of the torus is closer to the surface and the radial energy profile 
is flatter. It seems likely that here, magnetic diffusion alone cannot 
explain the evolution and that thermal diffusion is important. 

In an axisymmetric equilibrium, the buoyant-rising effect of 
thermal diffusion is resisted by the tension of the toroidal compo- 
nent. As iReisenegged d2007t) points out, this is likely to be more 
important when the torus is buried deep inside the star, which could 
explain why magnetic diffusion may be more important in deeply 
buried field 

As axisymmetric equilibrium with a flatter radial energy pro- 
file diffuses outwards, it experiences the same loss of toroidal flux 
as the non-axisymmetric equilibria. A purely poloidal axisymmet- 
ric field is unstable around the neutral line, i.e. where the poloidal 
field goes to zero. The instability, which drives poloidal field loops 
in a direction parallel to t he field's axis of symmetry , has recently 
been seen i n simulatio ns (Braithwaite & Spruit 2006, a nd in r otat- 
ing star s inlBraithwaitell2007l) . lMarkev & Tavled i ll 97311 19741) and 
Wright) d 19731) find that the instability can be suppressed by ad- 
dition of a toroidal field of at least ~ 1/4 of the strength of the 
poloidal component. This works because the instability increases 
the length of the neutral line, so it has to work against the tension of 
the toroidal component. The stronger the toroidal field, the smaller 
azimuthal wavenumbers are stabilised and just below the threshold, 
only the largest-scale (m = 2) azimuthal mode is unstable. This is 
because the larger wavenumbers stretch the toroidal field more than 
the smaller wavenumbers, so are easier to stabilise. Therefore as an 
initially stable mixed poloidal-toroidal field diffuses outwards and 
loses toroidal flux, it eventually becomes unstable to this m = 2 
mode. Until now, the non-linear development of this instability had 
not been studied. In light of the arguments in Section [3~2l it can 
be seen that once this threshold is crossed, it becomes energetically 
favourable for the flux tube to become longer and narrower. In other 
words, a falls from 180° to a smaller value, and the equilibrium 
becomes non-axisymmetric; in fact, exactly of the type described 
above. Since the analysis given in Section [3~2l assumed small a, 
it is not immediately obvious whether this transition takes place 
quasi-statically or whether a falls quickly from 180° to an equilib- 
rium at some lower value. We can appeal here to simulations, which 
appear to show a gradual transition on a diffusive time-scale. Fig- 
ure[TT]shows W, the length of the B r = line on the stellar surface 
(divided by the circumference) for an initially axisymmetric equi- 
librium which crosses this threshold due to toroidal flux loss. The 
rise in W does not seem sudden, but happens only on a diffusive 
time-scale. Figure[T"2lshows the field at three times: before, just af- 
ter, and some time after this transition, corresponding to //ta = 2.9, 
16.1 and 23.5 in Figure [TT1 

In addition, the instability of an axisymmetric poloidal field 
presents us with another way of understanding the adjustment, in 
a non-axisymmetric configuration, of a narrower flux tube to equi- 
librium described in Section [3~2l If the toroidal component in the 

4 This is equivalent to metal bar which is initially hot in the centre: heat is 
lost through the sides, but is also conducted along the bar so that at first the 
ends become warmer despite a fall in the total thermal energy. 

5 He also notes that in a star with a barotropic equation of state or which 
is otherwise non-stably stratified, a flux tube should rise to the surface on a 
dynamical time-scale. The issue of stable fields in non-stably-stratified stars 
will be looked at in a forthcoming paper. 
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Figure 11. The length of the B, = line on the stellar surface (divided by 
the circumference) for an initially axisymmetric equilibrium which crosses 
this threshold between axisymmetric and non-axisymmetric equilibria due 
to toroidal flux loss. 

flux tube is a little too weak, there is a long-wavelength instability 
around its axis, the resultant latitudinal displacements increasing 
the length of the flux tube and taking energy from, and weaken- 
ing, the poloidal component until stability is restored. Therefore as 
the toroidal flux is gradually lost into the atmosphere, the flux tube 
adjusts quasi-statically and is always marginally unstable. 



5 DISCUSSION 

In this section I consider some applications of these results to real 
stars. 

5.1 Comparison to observations 

As outlined in Section Q] both axisymmetric and more complex 
equilibria occur in all three classes of non-convective star: upper- 
MS, WDs and NSs. This implies either that the initial conditions 
vary - in particular the degree of central concentration of the initial 
field and/or flux connection through the star's surface, or that all 
stars form axisymmetric equilibria which then evolve after some 
time into more complex configurations via the process described in 
Section |4~4l Another way of looking at this is that there is a contin- 
uous series of equilibria and that diffusive processes move the star 
along this series. At the beginning of the series are deeply-buried 
axisymmetric fields, which evolve on the magnetic diffusion time- 
scale R 2 /n, where 77 is the magnetic diffusivity. As this happens, the 
field strength on the surface may increase at first and then decrease. 
After that are less deeply-buried axisymmetric fields with a greater 
fraction of their flux going through the star's surface; these evolve 
on a time-scale of either R 2 /n or R 2 /3/k where k is the thermal diffu- 
sivity; the field strength on the surface is likely to fall with time. Af- 
ter that come the non-axisymmetric fields which evolve on a time- 
scale of a 2 R 2 {3/K or a 2 R 2 /n; since a is continuously falling, so the 
time-scale of evolution falls until the field evaporates completely. 
During this phase, the field strength on the surface is likely to fall on 
a magnetic diffusion time-scale, but the dipole field strength could 
fall more quickly. The star begins its life somewhere on this se- 
quence and moves along with time. It is useful now to look at the 
actual time-scales in each of the classes of stably-stratified star. 
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Figure 12. The successive snapshots of an initially axisymmetric equilibrium which has evolved from that in Figure|2] A transition is made from axisymmetric 
to non-axisymmetric configuration as diffusive processes destroy toroidal flux and it becomes energetically favourable for the circular flux tube to lengthen 
and become narrower. The resulting equilibrium is equivalent to that reached directly from initial conditions with a flatter radial energy profile (Figure[3j- 



In an intermediate-mass MS star, the Alfven time-scale is 
~ 10 years for a field of lkG and the magnetic and thermal dif- 
fusion time-scales are FC-jr] ~ 10'° years and R 2 /k ~ 10 6 years; 
the lifetime of the star is typically 10 8 ~ 9 years. Composition gradi- 
ents in the radiative envelope are probably unimportant and since 
B ~ 10 6-12 , the buoyant rise time-scale - as given by \19\ - is 
greater than the lifetime of the star, unless the flux tubes are very 
narrow. Very narrow flux tubes would be found observationally 
by a very strong average field modulus compared to the averag e 
line-of-sight component, which is not seen (see e.g. lMathvsll200ll) . 
There should not therefore be any significant evolution of the field 
during the main sequence. Indeed, there is arguably no observa- 
tio nal evidence for fie ld evolution during the main-sequence (but 
sec iHubrig et alJboOOh . We need therefore different initial condi- 
tions. It is possible that stars are born with a range in central con- 
centrations around the threshold so that some stars develop axisym- 
metric fields and some more complex equilibria, but still with large 
a. The questions of the origin of the range in field strengths and the 
apparent cutoff at ~ 200 gauss (no stars are observed with fields 
weaker than this d own to an observational limit of a few gauss - see 
Aurie reet alj2007r) are still open, but since the former is generic to 
all non-convective stars and the latter is unique to A stars, it is likely 
that the solutions to these two questions are unconnected, and that 
there is some mechanism in A stars to destroy fields below a certain 
threshold. Whether there is a similar cutoff in O and early B stars 
as that in A and late B stars remains to be seen. 



In WDs, the Alfven, magnetic and thermal diffusion time- 
scales are a few days (field of 10 MG), ~ 10 10 and ~ 10 4 years 
respectively (Chanmugam & Gabriel 1972); Hall drift may also be 
important dMuslimov et al.l 1995). Because the thermal diffusivity 
is so high, the thermal buoyancy timescale given by d 1 9 lis less than 
the cooling time (10 9 ~ 10 yr) for the stronger fields. However, a com- 
position gradient can prevent buoyant rise because the species dif- 
fusion in WDs is very low. Since carbon and oxygen have almost 
exactly the same mass per electron, C-0 WDs lack composition 
stratification but the heavier O-Ne WDs are stratified. This could 
be the reason that we observe that th e magnetic WDs tend to be 
heavi er than the non-magnetic WDs (Wic kramasinghe & Ferrariq 
2005). Observationally, no correlations of magnetic properties with 
age have been found. As in A stars, a mix of dipole-dominated and 



quadrupolar and higher order fields are found, but none dominated 
by very high order harmonics. 

In NSs, the Alfven time-scale is ~ 100s for a field of 10 12 
G. The diffusive time-scales are subject to considerable uncertainly 
from both the theory - because Ohmic diffusion, ambipolar diffu- 
sion and Hall drift as well as superconductivity and superfluidity 
may all play roles - and from the observations. There is evidence 
that the fields of radio pulsars decay on time-scales of > 10 6 years 
while magnetar fi elds seem to decay faster, perhaps over around 
10 4 years (see e.g. Woo ds & T hompson 2004). Hall drift and buoy- 
ant rise provide natural explanations for this as they scale with the 
field strength, as B and B 2 respectively. Buoyancy of flux tubes in 
neutron stars depends mainly on the composition gradient, the dif- 
ference with WDs being that not only can chemical elements dif- 
fuse into or out of a flux tube but that the chemical elements can 
reach equilibrium via the neutrino-emitting reaction n <-» p + e 
jReisen egger 2007). The speed of the latter does not depend on the 
length scale, so fields with narrow flux tubes would evolve no faster 
than simpler fields if this were the dominant mechanism. Some 
NSs are observed with rather small-scale fields, such as IE 1207.4- 
5209 whose dipole field i s x 50 times smalle r than the typical fiel d 
strength on the surface jPavlov et all |2002| . ISanwal et all |2002|) . 
Such a field could be a low-o- equilibrium. This kind of equilibrium 
could provide an explanation for the observed phenomenon of NSs 
with similar dipole field strengths and rotation periods which have 
quite different observational propert ies, some being ob served as ra- 
dio pulsars and some as magnetars dMorris et al.f 2002). Magnetars 
with 10 14 G dipole fields could have a complex fields containing 
much more energy than the radio pulsars with the same dipole field 
strength, which could account for the high energy output of mag- 
netars. 



5.2 Permissible poloidal/toroidal ratios in axisymmetric 
fields 

In an axisymmetric configuration, it is only the poloidal component 
we can see on the surface and which manifests itself in various pro- 
cesses taking place outside the star, such as disc accretion, torques 
in a binary system, pulsar spindown, etc. and it is therefore only the 
poloidal component which we observe directly. The magnetic field 
inside the star, including the toroidal component, can be important 
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however in other ways, for instance as an energy source for mag- 
netars and as a means of creating 'hot spots' on neutron stars via 
anisotropic heat conduction. Also, the Lorentz force gives rise to 
a distortion of the density field which has two main consequences: 
firstly, by changing the moment of inertia and making the star tri- 
axial, free precession is expe cted. There is obse rvational evidence 
of precession in some stars jAkgiin et al1l2006h . Damping of this 
free precession results in a star with its rotation and magnetic axes 
either parallel or perpendicular, depending on whether the star is 
oblate or prolate, respectively. A spinning neutron star with non- 
par allel axes should be obse rvable via gravitational radiation (see 
e.g. iDall'Osso & Stellall2007l and references therein), especially a 
fast-spinning magnetar as required by 'millisecond magnetar' mod- 
els of SNe. An important question is therefore whether the field in- 
side the star is dominated by the poloidal or by the toroidal compo- 
nents, since the former makes the star oblate and the latter prolate, 
as well as because the decay of a large 'hidden' toroidal field can 
provide an energy source for magnetars. 

The energy argument of Section 13.21 as applied to axisym- 
metric fields gives us an important clue as to the permissible 
toroidal/poloidal ratios in a stable axisymmetric equilibrium: the 
poloidal field cannot be much stronger than the toroidal, or it 
would be energetically favourable to make a transition to non- 
axisymmetric equilibrium. However, it cannot be ruled out that the 
poloidal field is still strong enough that the star is oblate; more 
quantitative calculations are required. It also cannot be ruled out 
that the toroidal field can be much stronger than the poloidal; in this 
situation a non-axisymmetric equilibrium would adjust, the flux 
tubes becoming shorter and fatter, but an axisymmetric flux tube 
cannot become any shorter and fatter (a cannot go above 180°). 
The question of exactly how much stronger the toroidal fi eld can 
be ha s to do with an instability of the toroidal field itself l iTavlerl 
Il973h : this will be examined in a forthcoming paper. 



6 CONCLUSIONS 

Numerical magnetohydrodynamic simulations of a stably-stratified 
star have been performed with an initially random magnetic field. 
It is found that this initial magnetic field always evolves on the 
Alfven time-scale into a stable equilibrium configuration consist- 
ing of twisted flux tube(s). There are approximately axisymmetric 
equilibria with one flux tube forming a circle around the equator, 
and more complex, non-axisymmetric equilibria consisting of one 
or more flux tubes arranged in a more complex pattern with their 
axes lying at roughly constant depth under the surface of the star. 
Whether an axisymmetric or non-axisymmetric equilibrium forms 
depends on the radial profile of the initial field strength: a more 
centrally concentrated field evolves into an axisymmetric equilib- 
rium and a more spread-out field with greater flux connection to 
the atmosphere evolves into a more complex equilibriunj^. Fur- 
ther, higher-resolution simulations are required to quantify better 
this distinction in initial conditions, but it seems that in an ideal- 
gas star with polytropic index n = 3, if the initial field strength is 
tapered as B ~ p p then the threshold is p ~ 1/2. 



6 A magnetised region of initially 'turbulent' field, embedded in an ini- 
tially stationary homogenous conducting gas, is found to evolve into an 
axisymmetric equilibrium (this will be explored in a forthcoming paper.) 
It seems then that fiofi-axisymmetric equilibria form only because of the 
non-conducting nature of the atmosphere. 



These equilibria evolve quasi-statically due to diffusive pro- 
cesses (finite conductivity, thermal/composition diffusion). Ax- 
isymmetric equilibria rise towards the surface and eventually make 
a transition to the non-axisymmetric class; the non-axisymmetric 
configurations undergo a gradual lengthening and narrowing of 
their flux tubes. This is caused by loss of toroidal flux into the atmo- 
sphere. In configurations with more than one flux tube, each tube 
may have either positive or negative magnetic helicity; whether 
negative or positive has no effect on the stability. Stable zero- 
helicity equilibria are possible. 
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